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Abstract: The scaling and similarity laws concerning the propagation of isolated spherical blast 
waves are briefly reviewed. Both point source explosions and high pressure gas explosions are 
considered. Test data on blast overpressure from the interaction and coalescence of spherical blast 
waves emanating from explosives in the form of shaped charges of different strength placed in the 
vicinity of a solid propellant stack are presented. These data are discussed with regard to the 
scaling laws concerning the decay of blast overpressure. The results point out the possibility of 
detecting source explosions from far- field pressure measurements. 


1 Introduction 

The study of the propagation of blast waves in air is of great technical interest in many industrial 
applications- chemical, nuclear and aerospace. Blast waves are generated by point explosions 
(nuclear explosions and detonation of solid explosives, solid and liquid propellant rocket motors), 
high pressure gas containers (chemical explosions) and laser beam focusing. Blast wave problems 
also arise in astrophysics, hypersonic aerodynamics and hypervelocity impact. An understanding 
of the properties of the blast waves both in the near-field and the far-field is useful with regard to 
the characteristics such as blast overpressure, blast speed, and impulse. 

Many theoretical and experimental studies were reported by several investigators on the planar, 
cylindrical and spherical blast waves since the classical work of Taylor [1-3] on strong (intense) 
blast waves from point source explosions (Fig. la). Notable theoretical studies on point source 
explosions include those of Sedov [4-5], and von Neumann [6]. The similarity solutions were 
related to the description of strong shocks in the early stage of development covering short 
distances from the center of explosion. Blast waves of moderate shock strength emanating from 
point sources were theoretically studied by Sakurai [7, 8] with a perturbation scheme, Oshima [9] 
with a quasi-similarity method, and numerically by Goldstine and von Neumann [10], and Bach 
and Lee [11]. 

In the other extreme of the nearly linear regime, when the shock has become rather weak, 
Whitham [12] proposed a general theory to cover such flows. It describes attenuation of spherical 
shocks at large distances from the origin. The character of blast waves at large distances from the 
origin was also reported by Landau [13]. 

Numerical solutions to the more realistic case of initially isothermal high pressure gas 
explosions (Fig. lb) were provided by Brode [14, 15], Friedman [16], Liu et al [17], Chen and 
Liung [18]. The existence of secondary shocks under certain circumstances is demonstrated. For a 
detailed discussion on isolated spherical blast waves, see Sedov [5], Sachdev [19] and Kinney and 
Graham [20]. 

Interesting problems also arise in the propagation and interaction (coalescence) of blast waves 
emanating from the detonation of closely spaced explosive charges (sequentially detonating, 
spatially distributed ammunition stacks). Regions of significant overpressure associated with the 
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coalescence of blast waves from distributed ammunition stacks may be less extensive than 
regulatory requirements (Starkenberg and Benjamin [21]). 

This paper summarizes the characteristics and scaling laws associated with isolated spherical 
blast waves. Test data for the propagation and coalescence of blast waves from closely spaced 
multiple charges and solid propellant are examined and analyzed. 


2 Isolated Blast Wave Propagation 

A knowledge of the properties of isolated blast wave propagation is useful in the understanding of 
interaction and coalescence of blast waves. Here we briefly summarize the similarity nature of 
spherical blast waves originating from point-source (massless) explosions and high pressure gas 
containers. 

When disturbances of finite amplitude are propagated in perfect gases (destitute of viscosity or 
heat conductivity), discontinuities in pressure, velocity and temperature of the medium may occur. 
These are called shock waves (or shocks) or blast waves (intense shock waves). Shock waves and 
the flowfield become planar, cylindrically symmetric, or spherically symmetric, respectively if the 
energy source is in a plane, along a line, or at a point in space. The appearance of shocks is a 
consequence of nonlinear character of the equations governing the propagation of finite 
disturbances that is in part responsible for the formation of discontinuities. 


2.1 Boundary Conditions 


2.1.1 Rankine-Hugoniot Relations 


The Rankine-Hugoniot relations serve as boundary conditions at the shock front for the analytical 
or numerical solutions. Behind the shock, a strong blast wind is initiated, which moves in the 
direction of the shock. The Rankine-Hugoniot relations (boundary conditions) across a shock of 
finite strength at r = /?(r)are (Shapiro [22]; Liepmann & Roshko [23]): 


Pi = r + 1 
Po 



( 9 \ 

Pl 2 Y 

po r +1 Uo 


r- \ l h_ 

y + l’ U 


2 

y + 1 


/ 

1- 



U L 


(i) 


where p, p , and T denote respectively the pressure, density and temperature of the gas, y the 
specific heat ratio (isentropic exponent), cthe speed of sound, and u p the particle velocity 

immediately behind the shock. Here the shock velocity U (relative to the undisturbed gas; i.e. 
relative to a fixed coordinate system) is defined by 


U = dR/dt 


( 2 ) 


with R denoting the shock radius (or shock trajectory), and t the time. The speed of sound in the 
undisturbed medium is expressed by 


co = -Jm/po ( 3 ) 

The shock Mach number M s {-U I cq) and the particle velocity u (near the shock) may be 
expressed in terms of the pressure ratio as 
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The quantity P = (/?i / Po ~~ 0 represents the strength of the shock, and Ap = - p 0 denotes the 

overpressure. 


Fig. 2a shows the variation of shock Mach number and shock speed as a function of 
overpressure, as given by eq. (4a). The shock speed increases with increasing overpressure. For 
example, the overpressure corresponding to a shock Mach number of 3 is about 1 MPa. The 
corresponding shock speed is 1040 m/s, which is comparatively high relative to the ambient sound 
speed (335 m/s). The dependence of particle velocity on overpressure is provided in Fig. 2b, as 
obtained from eq. (4b). 


Strong Shock Conditions 

The strong shock conditions ( p x » p 0 , U » c 0 ), as obtained from ( 1 ), simplify to 
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In terms of the pressure ratio, the shock Mach number and the particle velocity may be expressed 
as 
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2.1 .2 Constancy of total Energy 

The energy released E a is assumed constant, and is expressed by 
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where the a = 0 , 1 and 2 respectively corresponds to the plane, cylindrical and spherical blast 
waves. This represents the explosion energy per unit area of the surface of the shock when 
R equals unity. In the strong shock limit ( Po « P )> eq. (7a) reduces to 


R R 

E = 4/rJ^-ipw 2 jr 2 dr + 4/rJ-“^-yr 2 dr (7b) 

0 0 

The constancy of energy in each sphere implies that the energy flowing into the sphere with the 
new material that enters it exactly balances the work which its original surface does by expanding 
it against the surrounding pressure. 
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2.2 Point Source Explosions 


2.2. 1 Strong Shock Solutions 

Self-similar solutions for the case of very strong shocks (short times after the explosion) were 
developed by Taylor [1-3] for spherical blast waves in his classical blast wave theory. These 
solutions are based on the form of the shock radius, which yield similarity solution. According to 
Taylor [2], a finite but large amount of energy is suddenly released by nuclear fission in an 
infinitely concentrated form. This leads to the formation of a shock wave according to 

4 )= s(rW ,5 E U5 t 215 = sir^E/pof 5 ( 8 ) 

which leads to a similarity solution. Here the quantity S is a dimensionless constant which is 
dependent on the specific heat ratio y , and is determined by solving the equations of motion. 

The governing partial differential equations (of continuity, momentum and energy) for a 
spherical wave in Eulerian form are thereby reduced to an ordinary differential equation (ODE), 
involving a single independent similarity variable defined by 

Z = r!R(t) (9) 

The Rankine-Hugoniot relations for strong shock (eq. 5) are utilized as the boundary condition at 
the shock front. Taylor [2-3] numerically solved the ODE. Taylor [2-3] solved the system for 
y - 1 .4 numerically, starting from the condition at the shock £ = 1 . 

An exact closed form analytical solution to the Taylor problem is given by Sedov [4], von 
Neumann [6], and later by J.L. Taylor [24]. Sedov [4] obtained solutions for plane, cylindrical and 
spherical blast waves, and for various values of y , Von Neumann [6] sought a similarity solution 
for the point explosion in Lagrangian formulation. The unknown shock boundary as a free surface 
is determined as part of the solution. 


2.2.2 Shocks of Moderate Strength 

As the shock strength decreases with time, there emerges a negative pressure phase (expansion 
region) behind the shock. Fig. 3 shows a typical pressure-time curve for the blast wave profile with 
a negative pressure region. The quantity tj represents the duration of the positive pressure phase. 
Exact numerical solutions (Sedov [4]) indicate an interesting feature connected with the 
development of a negative phase in the velocity distribution as the shock becomes weaker. 

The shock radius in the strong shock limit ( » P\ ) and the acoustic limit corresponding to 

short and large distances respectively are characterized by 

ns 

R oc t strong shock (short distances) (10a) 

R = c$t sound wave (at large distances) (10b) 


a. Sakurai ’s Solution 

For intermediate times when the shock strength becomes finite, small departures from the classical 
solution due to counterpressure effects are accounted for in the perturbation solution by Sakurai 
[7]. Sakurai [7] generalized Taylor’s similarity solution [2-3], and proposed a perturbation scheme 
for moderately strong shocks, in which the Taylor solution for strong shocks appears as the zeroth 
order solution (Sachdev [19]). Sakurai expressed the solution in terms of two variables (also see 
Sedov [5]], p. 297) 


4 = r/R(t\ T, = [c 0 /U{tf=ct/R 2 =UM1 (11) 
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The Rankine-Hugoniot equations for shocks of arbitrary strength (eq. 1) are utilized here as the 
boundary condition at the shock front. 


The similarity parameter of Sakurai [7] is expressed by 

c 0 t/R 0 =f{R/R 0 ) (12) 

where cq denotes the undisturbed speed of sound, R the shock radius, and Rq characteristic 
length. A quasi-similarity rule was proposed by Oshima [9]. 

b. Solution of Bach and Lee 

Bach and Lee [11] reported an approximate solution similar to that of Sakurai [25] for the entire 
range of the blast evolution (from strong shock to weak acoustic wave). The main assumption is 
that the density behind the shock wave in the blast obeys a power law, with the exponent being a 
function of time. A quasi-similarity is shown to exist for the decay coefficient 0 for moderate 
blast strengths as (Higashino et al. [26]) 

0 = 0(tj) (13) 

where 0s rr/r^A.^. j, n = MM] (14) 

The results of Bach and Lee [11] agree closely with the numerical solution of Goldstine and 
Von Neumann [10]. The condition 0 -- 1 corresponds to the strong shock limit, and 0 - Oto the 
acoustic (weak shock) limit. Sakurai’s perturbation solution (Sakurai [7-8]), and the approximate 
linear velocity profile solution of Sakurai [27] are shown to be less accurate. The velocity profiles 
suggest that the particle velocities can be negative in regions behind the blast wave in the moderate 
and weak shock regions (i.e., particles are directed towards the center of symmetry). These results 
are consistent with the exact numerical solutions of Sedov [5], 


2.2.3 The Shock Wave at Large Distance 

At large distances from the center of explosion when the shock has become rather weak, Whitham 
[12] proposed a general theory to describe the attenuation of spherical shocks. The only 
assumption is that since the entropy changes are of third order in shock strength, the flow with 
weak shocks may be considered isentropic. Thus R = /?(c 0 ,f)at large distances from the origin. A 
fraction of explosion energy which is degraded as heat is thus not available for doing work as the 
shock propagates. 


2.3 High Pressure Gas Explosions (Bursting Sphere) 

A more realistic problem than that of point explosion is the case of a unit sphere of high pressure 
gas at time / = 0 suddenly expanding into a homogeneous atmosphere (air) at rest. This problem 
was treated numerically by Brode [14-15] and analytically (approximately) by Friedman [16]. At 
sufficiently high pressure, an important new phenomenon arises here due to the occurrence of a 
secondary shock that does not appear in the one-dimensional shock tube problem (Sachdev [19]). 

A typical blast wave pressure profile is shown here (Fig 4). The major difference between a 
planar shock wave and a planar blast wave is that the flow properties behind the shock wave front 
are constant. A distinct characteristic of a spherical blast wave is a moving shock wave (jump in 
pressure) immediately accompanied by expansion waves behind the blast wave front. There is a 
blast wave profile region in which the flow properties are rapidly decreased due to expansion 
waves (Chen & Liang [18]). Due to its 3-D expansion, the intensity of the blast wave always 
decreases with time and distance during its propagation. The overpressure phase is due to the 
shock front and the under-pressure phase is due to the expanding waves following the shock front. 
Both phases result in two visible shock peaks on the pressure-time plot. Because the blast wave 
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intensity decreases with time, the range ( r t < r < /y ) of the blast wave profile increases, and the 

magnitude of the pressure peaks (maximum and minimum) of the blast wave decreases (Chiang 
[18]). By contrast in shock tube experiments, the pressure downstream of a shock wave is often 
constant between itself and the following contact surface. 

2.3.1 Brode’s Numerical Solutions 

Brode [14] simulated spherical blast waves, and considered two cases of practical importance: 1. 
strong shock, point-source explosions, 2. the sudden release of high pressure isothermal (both cold 
and hot) gas from spherical enclosures. For the point explosion problem, he considered the exact 
solution of von Neumann [6] as an initial condition. Ideal gas equation of state is considered. 
Artificial viscosity, introduced by von Neumann and Richtmyer [28] is considered for stabilizing 
the numerical solution. 

Brode [15], in a subsequent study, considered a blast wave from a spherical charge of 
trichlorotoluene (TNT). An equation of state applicable for TNT is utilized. In his numerical 
studies of detonation with a spherical charge of TNT, the existence of a second shock in the focal 
region was demonstrated. By considering real gas at high temperature, Brode [15] has shown that 
the blast over pressure at any radius is lower than that given by the ideal gas solution, since a part 
of the available energy is absorbed by the processes of ionization and dissociation. 

Brode [14] and Goldstine and von Neumann [10] show that the variation of shock overpressure 
varies with the shock radius R follows a power law of the form 

p-\ = AR n (15) 

where n is a slowly varying function of R , and A a constant (Sachdev [19], p. 242), where p is in 
atmospheres. 

Brode [15] has shown theoretically that the peak overpressure p sQ from a free spherical air 
blast can be expressed in dimensionless form (Wyle Labs [29]) 
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where p so is the peak overpressure due to TNT explosion, p 0 the atmospheric pressure, A the 
reduced distance, r the actual distance, and E the energy of the blast. 


2.3.2 Numerical Solutions of Liang, Wang and Chen 

Recently Liang et al. [30] numerically investigated the character (propagation and decay) of a 
spherical blast wave emanating from an isothermal bursting (rupture of a) sphere of pressurized air 
in a free field. Three minor shock waves moving behind the main shock are observed. 


2.4 Variable Energy Blast Waves and Blast Wave Anisotropy 


2.4. 1 Variable Energy Blast Waves 


Variable energy blast waves present themselves in spark discharges, laser induced blast waves and 
two-phase detonations. The instantaneous constant energy blast wave represents a special case of a 
more general class of variable energy blast waves in which the energy input varies with time 
according to a power law of the form [31] 


E = E 0 t? 
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where E 0 and p are positive constants. According to similarity theory, the shock position under 
these circumstances varies as 


R~ft- + P)'S (17b) 

The case of P - 0 reduces to the instantaneous blast wave, whereas the case of p * 0 corresponds 
to the expanding piston problem. 

Freeman [31] investigated in detail the special case of linear energy input (p = 1 ). The data of 
Freeman [31] suggest that the shock wave radius at early times is smaller than that according to 
instantaneous energy blast wave theory [32]. 

2.4.2 Anisotropy in Blast Waves 

Merlen and Dyment [33, 34] investigated similarity solutions for the anisotropic blast waves in a 
stationary ambient, and explosions in a moving gas with applications to gun-firing and the 
associated muzzle wave propagation. The anisotropy can arise on account of initial conditions 
(energy release or distribution) or external flowfield (moving stream) [33, 34]. In a stationary 
ambient, the anisotropy is occasioned by non-instantaneous release of energy (variable energy) or 
the azimuthal distribution of energy. Anisotropy can also occur with an isotropic explosion in a 
moving stream. 

a. Intense Anisotropic Blast Waves with Azimuthal Energy Release 
An axisymmetrical azimuthal distribution of energy may be represented by [34] 

E = E\e{(f)t & §<<j)<7t (18a) 

where ^ denotes the azimuthal angle. If the energy release is not instantaneous (/? * 0 ), a vacuum 
forms inside the flow, and the empty volume will be aspherical [34], The case of 
p = 1 approximates the fire-arm shot (gun firing). The shock radius is expressed in the similarity 
form 
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It was shown that [34] in the farfield the wave properties depend only weakly on the anisotropy of 
the energy release, and depends primarily on the total energy, so that one can regard the wave as 
spherical with an equivalent energy, 

b. Isotropic Explosion in a Moving Gas 

The effect of ambient pressure /? 0 (counter pressure effect) and ambient velocity (t/ 0 ) have been 
investigated by Merlen and Dyment [32] with the aid of a perturbation scheme relative to the self- 
similar solutions of isotropic point explosions. The anisotropic effect of ambient motion and the 
isotropic effect of ambient pressure are examined. The solutions were sought in the form 
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and is a constant. 
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The theoretical solutions were compared to the measurements on gun firing with various gun 
calibers [33, 34]. It was demonstrated that the axial propagation of the blast waves X emanating 
from a gun can be expressed in dimensionless form [33] as 

j~Po~^Q~ X 

ypg c s D 

where the subscripts 0 and g respectively denote the ambient gas conditions and the gas at the gun 
muzzle. The quantity D denotes the gun caliber. 

The muzzle wave tends asymptotically towards the blast wave of a non-instantaneous intense 
point explosion [33]. It was shown that the velocity of the external flow reinforces the shock 
towards the upstream and weakens it downstream. The pressure of the external medium becomes 
important only after the shock weakens, as is to be expected. 

A generalization of the muzzle wave similarity was later reported by Merlen [35]. This work 
describes an improvement of the muzzle wave similarity that permits treatment to unsteady and 
non-sonic gun conditions. A numerical study of blast flowfields induced by supersonic projectiles 
discharged from shock tubes was published by Jiang et al. [36]. 
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3 Interaction of Spherical Blast Waves (Blast Wave 
Coalescence) 

Reported studies on interacting spherical blast waves are relatively few. The interaction of a pair 
of (weak) cylindrical blast waves of either equal or unequal strength were studied experimentally 
and computationally by Higashino et al. [26]. Liang et al. [30] numerically investigated spherical 
blast wave propagation and reflection. In sequential detonation and spatially distributed blast 
sources, it is known that shock waves coalesce and focusing can lead to significant amplification 
of peak overpressure. At sufficiently large distances from the charges, coalescence will occur. An 
overtaking wave may first encounter the negative overpressure phase of the wave being overtaken, 
causing it to decelerate. The overtaking wave must be somewhat stronger than the wave being 
overtaken or it will not be able to penetrate the negative overpressure phase and coalescence 
cannot occur. Once this phase has been penetrated and the overtaking wave has penetrated the 
positive phase of the wave being overtaken, coalescence takes place but occurs only after some 
additional propagation. 


3.1 Peak Overpressure 


Several approximate methods of combining peak overpressure from N coalesced waves are 
described in Starkenberg and Benjamin [21] as follows: 


2.1.1 Simple superposition 

N 

a p = ( 21 ) 

i=i 


where A p represents the combined overpressure, and A p t the overpressure of blast wave / . 

3.1.2 Full LAMB Model 

Hikida and Needham [37] utilized the following scheme for blast coalescence: 
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where p refers to density, and U the velocity of blast propagation. 


3.1.3 Single Effective Charge 

The method (conception) of single effective charge considers the use of pressure produced by the 
total explosive weight of all stacks with the coalesced waves combined at their centers of charge. 
That is, the peak overpressures associated with coalesced waves are essentially the same as those 
produced by single charges of the same total explosive weight. Zaker [38] experimentally showed 
that with three equal charges, the third pulse tends to overtake the second before the second 
overtakes the first. Thus 


E 


tot 



(23) 


Zaker [38] experimentally showed that with three equal charges, the third pulse tends to overtake 
the second before the second overtakes the first. 

Charges of equal/unequal mass (charges) and delay time, and separation distances need to be 
considered in the development of pressure combination algorithms and coalescence maps 
(patterns). 


3.2 Structure of Wave Interaction 

According to the classification by von Neumann [39], the interaction of shock waves is called 
regular reflection when they first encounter each other (see Higashino et al. [26]). Later on, the 
interaction changes to an irregular type of reflection called by von Neumann as Mach reflection . 
Colella and Henderson [40] have shown that for very weak shock waves, the von Neumann theory 
yields no physically acceptable solution (the so-called von Neumann paradox ), and this type of 
weak irregular reflection is termed a von Neumann reflection. 


4 Test Data on Overpressure 

Blast overpressure data were obtained from scale model and full scale solid rocket propellant 
motors (upper stage) with embedded charges at NASA White Sands Test facility (WSTF) [41]. 
The solid propellant is struck by shaped charges of 500 gm and 20 gm. Both symmetrical (equal 
strength) and unsymmetrical (unequal strength) wave interactions were studied. Additional 
discussion of the test data is presented in [42]. 
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4.1 Scale Model Tests 


4.1.1 Test Configuration and Instrumentation 

In one of the scale-model tests, considered here, the propellant weight was 215 kg. For 
comparison, the propellant in the full scale motor is about 2,007 kg. The propellant composition is: 
71 % ammonium perchlorate, 18 % aluminum powder, 11 % hydroxyl-termicrated polybutadiene 
binder (HTPB). The configuration (Fig. 5a) includes a Ti-cylinder around the propellant stack to 
simulate the effects of the casing around the actual motor. The solid propellant was impacted with 
one 500-gm conical-shaped charge (CSC) of C4 and 17 gm of Comp A-4 booster pellet. The blast 
overpressure data is recorded up to 45 sec. from the fire command ( t = 0 sec). Only the near-field 
blast pressure was measured. 

The blast overpressure was measured by piezoelectric transducer pressure gauges (PCB) 
mounted on three adjustable elevated gauges lines (EGL), designated F, G and H, separated 120 
deg. Apart (Fig 5b). Five pairs of PCB’s were flush mounted equal distances apart on each side of 
the pipe, in designated positions A and B. The EGL’s were constructed from 101.6 mm diameter 
steel pipe, 9.144 m long, and were held 1.83 m above the ground (parallel to). The PCB’s had a 
500 kHz resonant frequency, and a 1- ps rise time. The PCB’s were located 180 deg from each 
other, parallel to the ground. EGL-F and EGL-H are fitted with a flat plate at one end, with the 
PCB mounted in the plate center, thus measuring the stagnation overpressure, which includes the 
dynamic components of overpressure due to the particle motion (impact). 

In compressible flow the isentropic stagnation pressure p t always exceeds the sum of the 
local static pressure p and dynamic pressure q , and is expressed by 


where 



(24a) 


(24b) 


q = ^pu 2 = ^yM 2 p 


(24c) 


Eq. (24a) suggests that in incompressible flow ( M 0 ) the dynamic pressure represents 
exactly the impact effects, i.e., the difference between the stagnation pressure and the static 
pressure (Kinney [20], Shapiro [22]). In compressible flow, the stagnation pressure is not simply 
the sum of static and dynamic pressure as in the case of incompressible flow. 


4.1.2 Test Data for Overpressure 

Fig.6a shows the overpressure history for the stagnation overpressure from the closest gauge on 
EGL-F (channel-67), which is 1.50 m from ground zero. This typical blast waveform is 
characterized by a peak overpressure of 0.207 MPa at 16.2 ms, and positive pressure duration 
(about 1 ms). We also notice a distinct wave due to propellant explosions, with a peak 
overpressure of 0.113 MPa at 19.2 ms. Only one major peak overpressure due to the propellant 
bum is observed here. The peak stagnation overpressure due to the charge is higher than that due 
to that of the propellant. Calculations (eq. 24) show that the estimated peak static overpressure at 
this location due to the charge is about 0.136 MPa. It is notable that the propellant pulse is 
somewhat wider than the pulse resulting from the 500-gm CSC. The relatively longer decay of the 
overpressure from the propellant blast is attributable to the long duration over which the energy is 
added to the flowfield. The overpressure phase is associated with the shock front location, and the 
under-pressure phase is connected with the expansion waves following the front. The data also 
suggest that the interaction between the blast waves from the 500-gm charge and the propellant is 
not appreciable. 
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The overpressure history for a number of locations along the EGL: F-A is demonstrated in 
Fig. 6b. This plot illustrates the propagation and attenuation of the blast wave. The peak 
Ap diminishes from 0.207 MPa at 1 .494 m to 0.0172 MPa at 7 m. Considering the first two PCB’s 
(separated by about 0.61 m), the 500-gm peak is displaced by about 1 ms. This result is consistent 
with the shock speed of about 457.2 m/s (Fig. 2a) at Ap « 0.1 1 MPa of the second peak. 


4.1.3 Peak Overpressure Extrapolation to the Far-Field 

Fig. 6c presents an extrapolation of the near-field peak A p data to the far- field. The near- field data 
are correlated by the power law relation 


Ap ~r 


- 1.06 


This result compares to the acoustic limit expressed by 

Ap ~ r -1 


(25) 


(26) 


4.2 Full Scale Propellant Test 


4.2. 1 Test Configuration and Instrumentation 

The full scale tests are described in Ref 41. Fig. 7a shows a photographic view of the actual 
motor. The overall dimensions of the test article were 1 .524 m wide x 2.286 m tall with 2007 kg 
of propellant. The propellant composition was: 71 percent ammonium perchlorate, 18 percent 
aluminum powder, 11 percent hydroxyl-termicrated polybutadiene binder (HTPB). The motor 
casing was fabricated from titanium alloy TI-6A1-4V. A total of four conical shaped charges 
(CSC) were used. An aft enclosure was used to simulate the nozzle. Two 500-gm shaped charges 
used for the breakup system were oriented 180° from each other. Two 20-gm shaped charges were 
located in the spin table assembly, near the top of the nozzle skirt. The 20-gm charges were 
located at azimuths of 191.3° and 338°, while the 500-gm charges were at 90° and 270°. The 500- 
gm CSCs are initiated 0.87 ms after the 20-gm shaped charges are actuated. 

The test set up and instrumentation was essentially similar to those described for the scaled 
propellant tests (Fig. 7b). The instrumentation for the full scale testing was similar to that in the 
scaled tests, except the charges were not symmetric for EGL-F and EGL-G. 


4.2.2 Test Data for Overpressure 

Fig. 8a illustrates the temporal history of the stagnation overpressure on EGL-F (channel-67), 
indicating a series of blast waves. Two 20-gm charges are followed by two 500-gm charge peaks 
in succession. The peak Ap of about 0.448 MPa due to 500-gm is considerably higher than that 
due to 20-gm peaks, since the energies are different, as indicated by the scaling law. The estimated 
peal static overpressure (from eq. 24) is about 0.348 MPa. For the EGL-F the two 500-gm are 
separated by about 1.22 m, and travel at about 701 m/s, so that At » 1.7 ms (the measured A t is 
about 1.7 ms). The peak Ap related to propellant explosion is superposed on to the 500-gm charge 
waveforms, with relative peak magnitude of about 0.138 MPa. This value of the propellant 
overpressure of 0.13 MPa for the full scale motor is of the same of order of magnitude as that in 
the scaled test. This suggests that the propellant fragment in both the scaled or full scale test is 
perhaps of comparable size. The existence of multiple peaks in the overpressure history represents 
the interaction between the blast waves emanating from the 500-gm charge and the propellant. 
Also secondary shock waves are noted. The 20-gm charges do not seem to have initiated any 
significant propellant bum (or breakup), with more localized rather than large-scale action. 

It is remarkable that the positive pressure duration of the peak of the first 500 gm charge is 
considerably less than that due to the peak of the second 500 gm charge, even though the peaks for 
the two charges are essentially the same. 
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The peak A p history for various locations along the EGL-F (Fig. 8b) suggests that the 500- 
gm charges ultimately coalesce at about 20 sec. This merging of shocks is a consequence of the 
fact that the trailing shock travels relative to the particle motion occasioned in the wake of the 
leading shock; thus it moves at a higher speed for the same strength. The A p decreases from 0.448 
MPa to 0.052 MPa as the distance increases (due to spherical spreading). The strength of the 
merged shock is due to the two charges. There seems to be a problem (data drop out) with channel 
37 after about 12 sec. 

The characteristics of the blast propagation for EGL: F-B transducers, shown in Fig. 8c, are 
seen to be similar to that of EGL: F-A in all respects. This result is to be expected, as the 
transducers are similarly situated, and demonstrates the consistency of measurements. Likewise, in 
the case of EGL:G-A (Fig. 8d), the behavior is similar to that of EGL-F, since the charges are 
similarly situated. The relatively small positive pressure duration of the peak of the second 500 gm 
charge to that of the first 500 gm charge is also manifested in the data of F-B and F-G also. In the 
case of EGL:H-A (Fig. 8e), the two 500-gm charges are located equidistantly relative to the 
gauges, so that the two peaks are spaced very closely. However they are located farther from the 
stagnation sensor (ch-67). Thus even though coalescence might have occurred, the measured 
stagnation pressure at this location is seen to be lower than that seen for EGL:F-A. 

4.2.3 Spectral Distribution of Energy 

The spectral distribution of energy for the pressure signal from channel 67 (F-A) for individual 
explosions of short duration for the various explosions (20 gm charge, 500 gm charge, and 
propellant explosion) are depicted in Fig. 8f. The quantity SPL denotes the sound pressure level 
(dB) relative to a reference pressure p re j of 20 //Pa , where 

SPL = \og w {plp ref } (27) 

It is seen that the energy is distributed over a broad range of frequencies, as is to be expected for 
impulse type functions typical of the blast wave. There is some spectral similarity for the various 
explosions, with the spectra displaying a distinct peak. The spectra are similar to those presented 
for blast waveforms (Diaci and Mozina [43]). 

4.2.4 Peak Overpressure Extrapolation to the Far-Field 

Fig. 8g presents the extrapolation of peak static overpressures to the far-field. The near-field data 
could be correlated by the power law relation 


Ap ~ r -0 ' 854 (28) 

The extrapolation suggests that the peak pressures at about 7.62 km are about 0.001 MPa, and thus 
are discernible. This consideration points to the feasibility of a ground microphone system for the 
possible detection of unanticipated (unexpected) onboard explosions. However, the occurrence of 
multiple explosions, blast wave coalescence, and absorption of sound at high frequency and over 
long distances should be taken into account, and present considerable complexity in the detection 
process. 


4.3 Scaling of the Farfield Overpressure 

The anisotropy of the blast waves represents an important consideration linked to farfield detection 
of explosions. The present data (pressure wave forms) suggest that in general the explosions are 
closer to instantaneous blast waves rather than variable energy blast waves. In the absence of a 
moving medium surrounding the initial blast, the anisotropy effects in the present experiments thus 
likely arise primarily due to the distributed energy sources (initial energies). As the wave 
properties in the farfield depend primarily on the energy levels of explosions rather than their 
anisotropic variations [34], it seems reasonable to regard the wave as essentially spherical with an 
equivalent energy. 
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With regard to the possibility of detecting source explosions from farfield measurements, the 
data concerning the spatial variation of peak overpressure from both the scaled model test and the 
full scale test (as shown in Figs. 6c and 8g) are reformulated in dimensionless form in an attempt 
to determine the existence of any similarity behavior. 


Fig. 9 depicts a plot of the overpressure vs. distance in dimensional form (corresponding to 
point source explosion) as 


A P 


=/ 


Po 1 _{E/p 0 ) ] 


1/3 


(29) 


In plotting this data, the total equivalent TNT of the propellant and the explosives have been 
considered. In particular, the equivalent TNT for the solid propellant is taken as 10 percent by 
mass, and for the 500 gm charge an equivalent of 0.707 kg of TNT is considered. 

The results suggest that both the scale model data and the full scale data for the peak 
overpressure match satisfactorily, considering a scale ratio of about 10 for the propellant mass in 
the present data. It is thus seen that the anisotropic effects on the farfield are not appreciable in the 
present experiments. This result tends to point out the feasibility of detection of source explosions 
from farfield measurements. Further experimental data are needed to confirm this assertion. 


5 Conclusions 

Both the scaled test data and full scale propellant blast test data have been helpful in our 
understanding of the blast propagation from single and multiple charges. The data illustrate the 
effects of blast coalescence and propagation of blast waves of different strength. The test data for 
blast overpressure are found to be generally consistent with the available blast scaling laws. In this 
test data, the coalescence effects do not appear to be significant with regard to the blast 
overpressure. Also the peak over pressure seems to follow a trend of nearly linear decay with 
distance (acoustic type) at relatively short distances from the explosion source. The scale model 
and the full scale data for the peak overpressure demonstrate a similarity nature. 
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point source shock front 



ambient gas at rest 
(po > Po * T 0 ) 


Fig. la Schematic of a point source explosion 


high pressure gas shock front 



Fig. lb Schematic of an explosion of an initially spherical charge. 
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Fig. 2a Shock velocity and shock Mach number from Rankin-Hugoniot relations. 
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Fig. 2b Particle velocity for blast wave propagation. 
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Fig. 3 Pressure-time variation of a blast wave of moderate shock strength 
(adapted from Kinney and Graham [20]. 



Fig. 4 A typical blast wave from a spherical charge, showing secondary shock. 
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Fig. 5a Test set up for scaled testing. 



Fig. 5b EGL and PCB locations for scaled test. 
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Fig. 6a Overpressure history at the stagnation location for the scaled test. 
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Fig. 6b Overpressure history for the scaled test (contd.). 
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Fig. 6c Peak overpressure scaling for the scaled test data. 
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Fig. 7b EGL and PCB locations for full scale testing. 
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Fig. 8a Overpressure history at the stagnation point for full scale testing. 
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Fig. 8b Overpressure history of sensors along F-A for full scale testing. 
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Fig. 8c Overpressure history of sensors along F-Bfor full scale testing. 
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Fig. 8d Overpressure history of sensors along G-A for full scale testing. 
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Fig. 8e Overpressure history of sensors along H-A for full scale testing. 



Fig. 8f Spectral distribution for full scale testing. 
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Fig. 8g Peak overpressure scaling for full scale testing. 



Fig. 9 Dimensionless peak overpressure variation. 
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